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ABSTRACT 



The number of rich galaxy clusters per unit volume is a strong function of O, the 
cosmological density parameter, and erg, the linear extrapolation to z = of the density 
contrast in 8/i _1 Mpc spheres. The CNOC cluster redshift survey provides a sample 
of clusters whose average mass profiles are accurately known, which enables a secure 
association between cluster numbers and the filtered density perturbation spectrum. 
We select from the CNOC cluster survey those EMSS clusters with bolometric 
L x > 10 45 ergs -1 and a velocity dispersion exceeding 800 km s" 1 in the redshift ranges 
0.18-0.35 and 0.35-0.55. We compare the number density of these subsamples with 
similar samples at both high and low redshift. Using the Press-Schechter formalism and 
CDM style structure models, the density data are best described with ag — 0.75 ± 0.1 
and f2 ~ 0.4 ±0.2 (90% confidence). The cluster dynamical analysis gives O = 0.2 ±0.1 
for which erg = 0.95 ±0.1 (90% confidence). The predicted cluster density evolution in 
an = 1 CDM model exceeds that observed by more than an order of magnitude. 



Subject headings: galaxies: clusters, cosmology: large-scale structure of universe 



1. Introduction 



The clustering of galaxies grows via gravity from density perturbations which are characterized 
by their power spectrum, P(k). Various theories predict the shape of P(k), but they do not 
accurately predict its amplitude. An integral constraint on the normalization is conventionally 
parameterized as as, the fractional mass variance in 8 hr 1 Mpc spheres calculated using the 
linear extrapolation of P(k). Rich galaxy clusters are particularly sensitive probes of as- N-body 
simulations have established that the Press-Schechter formalism (Press & Schechter 1974, hereafter 
PS) gives a remarkably accurate prediction of the number of clusters per unit cosmological volume 
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as a function of mass. Modeling the low redshift data with the Press-Schechter formula, and using 
CDM style P(k), leads to a range of possibilities, from = 1 and as — 0.5 to ~ 0.2 and as — 1 



and values that interpolate between the two (Henry & Arnaud 1991, White, Efstathiou & Frenk 



1993] [Eke, Cole fc Frenk 1996| , [Viana fc Liddle 1996j , Pond fc Myers 1996| ). The O dependence can 



be disentangled from as if data giving n(M) dM are available as a function of redshift ( Pukbir"fc 
Blanchard 1996|) . 

The Canadian Network for Observational Cosmology (CNOC) cluster sample and observational 
strategy ( |Yee, Ellingson Sz Carlberg 1996|) was specifically designed to produce data useful for a 
as measurement. The sample's primary advantage is that the cluster masses are accurately known 
near the virial radius which is essential for a reliable estimate of the linear mass scale from which 
the cluster collapsed. Here we combine our results with similarly selected clusters at higher and 
lower redshifts in Section3. The data are modeled in Section 4 to draw conclusions about the 
values of as and 0. 



2. Press-Schechter Predictions 

The number density of clusters in the mass range M to M + dM is predicted to be 

-(M)dM = d ^ exp Hg(«)/2A»] dM, (1) 

where 5 c (z) = 5o(Q)/D(z,Q) gives the linear overdensity at which a collapsed structure is 
approximately virialized. The function 5q = 0.15(127r) 2/,3 O a0185 is nearly constant at 5q ~ 1.68 
( Navarro, Frenk fc White 1996| ). The growth factor, D(z,£i), gives the redshift dependence 



of the linear amplitude of the density perturbations ( [Peebles 1993 ). The quantity A(r^) 



measures the fractional linear mass variance in spheres of radius rj,. It is calculated using a 
tophat filter from a parameterized version of the CDM spectrum (r fixed at 0.2 [Efstathiou 



Bond fc White 1992] ) whose normalization is adjusted such that as = A(8/i x Mpc). The 



"just virialized" radius is near l.Bh 1 Mpc and has an overdensity of ~ 1780 0,6 ( [White, 



Efstathiou & Frenk 1993; ). Near 1.5/j" 1 Mpc M(r) oc r p , where p ~ 0.64 for the CNOC 



data flCarlberg, Yee fc Ellingson 1997| ). The tophat filtering scale, rz, is related to Mi. 5, as 
r L ~ 8.430(^) a2p /( 3 - p )[Mi. 5 /6.97 x lO 14 ^" 1 M /O] 1 /( 3 " p )(l + z)-^ 3 ^/^ 1 Mpc (a generalization 
of [White, Efstathiou fe Frenk 1993| ) where fi(z) is the value of at redshift z. 



We integrate Eq. ^| from the minimum mass in the sample to infinity to derive n(> M1.5). 
For comparison with measurements, we average the density predictions over the redshift ranges of 
interest. These predictions are plotted in Figures [l| and [2] for O = 0.2 and = 1, respectively. 
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3. Cluster Density Estimates 



The CNOC sample was drawn from the EMSS cluster survey ( Henry et al. 1992| ) 



supplemented with a few high redshift clusters identified later flLuppino fc Gioia 1995| ). We 
impose the constraints f x > 4 x 10~ 13 erg cm -2 s^ 1 , L z (0.3 — 3.5Kev) > 4 x 10 44 ergs _1 , redshifts 
between 0.18 and 0.55, and in the declination range —15 to +65 degrees. We have recalculated the 
f x and L x to better allow for the redshift dependence of the flux in the "detect cell" . The main 
difference is that the cluster MS0302+16 now falls below our luminosity cut, although this has a 
negligible affect on the results (see also Nichol, et al. 1997| ). 



The Mi. 5 masses for each cluster in the CNOC sample are estimated to be 
M v bMv (l-5/i -1 Mpc/r v ) p where bj^v is the virial mass bias, r v is the virial radius, calculated as 
a ringwise potential flCarlberg et al. 1996| ), and M v is the resulting virial mass. For the CNOC 
clusters buv — 0.85 on the average ( |Carlberg, Yee fc Ellingson 1997| ), that is, the virial mass is 
always an overestimate. The masses Mi. 5 are calculated from the M v and r v of the CNOC clusters 



(Carlberg et al. 1996 as updated slightly by the use of the finalized catalogues). 



3.1. Correcting for X-Ray Selection 

We need to count the number of clusters having Mi. 5 larger than a specified mass within some 
redshift range. A significant complication is that we have a sample defined by its X-ray properties, 
but we want to make a measurement based on its distribution of characteristic masses. It is well 
known that there is a strong L x — Mi. 5 correlation, but with a substantial scatter (for instance, 



Edge Sz Stewart 1991 , hereafter ES). That is, some clusters that are in the L x limited sample will 
be below the specified mass limit and vice versa. To correct from an X-ray selected sample to a 
mass selected sample we proceed as follows. We use the Edge and Stewart study of the X-ray and 
optical properties of nearby clusters, which appear to be essentially identical to the CNOC sample 
in the relevant parameters, see Figure |3[ The ES sample is effectively a mass selected sample, at 
least for the rich clusters which are our concern here. All the selected clusters have a bolometric 
L x > 10 45 ergs -1 . Examining Figure [3|, we see that setting the minimum mass equal to that for 
a cluster with a v ~ 800 km s -1 will help maximize the number of clusters in the sample while 
keeping the corrections from the L x selected sample to the a v selected sample relatively small. 

There are two different correction factors for L x defined samples, depending on whether or 
not the parent sample has known velocity dispersions. In the ES sample there are 12 clusters with 
a v > 800 km s -1 , of which 6 are also above L x > 10 45 ergs" 1 . Therefore we estimate that for a 
sample selected with both the L x and a v limits, the true density of clusters with a v > 800 km s" 1 
is a factor of f xa = 12/6 = 2 ± 1 times the density measured. There are 8 ES clusters with 
L x > 10 45 ergs _1 for all velocity dispersions and 12 with a v > 800 km s -1 . Therefore for a sample 
selected with only the L x limit, we estimate that the number of clusters above a v = 800kms _1 
is fx = 12/8 = 1.5 ± 0.7 times the measured density. Although these corrections are not very 
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accurate, they illustrate that the corrections are not large, and their errors are comparable in size 
to those from the subsamples themselves. 



3.2. The CNOC Sample 

To convert from EMSS 0.3-3.5Kev luminosities to the bolometric luminosity we derive a 
mean bolometric correction of a factor of 2.20 at 6.8 Kev and 3.12 at 13.6 Kev, with only a small 
dependence on the HI column. We adopt a uniform bolometric correction of a factor of 2.5 for 



the CNOC clusters (see the high L x sample of Henry, Jiao & Gioia 1994 for representative EMSS 



cluster temperatures). Imposing the limit a v > 800 km s 1 reduces the CNOC sample size from 12 
clusters to 8 (2 of the 16 clusters observed being dropped as being below the X-ray limits, 1 cluster 



not from EMSS, and 1 has a poorly determined a v ). Using the V e /V max method (|Avni fc Bahcall 



1980 ) we then measure the volume density of these clusters, allowing for the EMSS sky area in the 



CNOC region as a function of flux (see Henry et al. 1992 for more details on this procedure and 



also an example table of EMSS all-sky coverage). The mean densities of these clusters for Q = 0.2 
are reported in Table |j. The CNOC sample is split into a low redshift subsample, 0.18 < z < 0.35, 
which has a smallest Mi. 5 of 4.8 x lO 1 ^ -1 Mq and a moderate redshift subsample, 0.35 < z < 0.55, 
which is found to have a smallest Mi. 5 of 6.7 x 10 14 /i _1 Mq. The f xa corrected densities are given 
in Table |l|. 



3.3. Low-Redshift Samples 

At low redshift there are several samples to consider as sources of density estimates. The 
most straightforward dataset is that of Henry & Arnaud (1991, hereafter HA) who provide an 
X-ray luminosity function at low redshift, which when integrated from L x = 10 45 ergs -1 to infinity 
yields an (f x corrected) volume density of ~ 7.5 x 10~ 7 /i 3 Mpc -3 . The velocity dispersions of these 
clusters extend below 800 kms , but none below 750 kms^ flZabludoff, Huchra, fc Geller 1990|) , 



which we adopt as the minimum velocity dispersion. We scale the cluster masses with velocity 
dispersion as M oc cr 3 with a reference value of 5.7 x 10 14 /i _1 M . We estimate the Mi. 5 of the HA 
sample as 4.7 x 10 14 /i _1 M Q . 



The ESO Cluster Survey ( Mazure et al. 1996 ) finds a cluster density of 2.5 x 10 h Mpc 



for a v > 800 km s -1 at z < 0.1. However, the ESO and CNOC velocity dispersions are not 
calculated in the same manner. The CNOC velocity dispersions are estimated using an explicit 
background subtraction. On the average we find that our velocity dispersions are about 7% lower 
than those calculated from precisely the same velocities using the iterated bi-weight estimator 
( Beers, Flynn fc Gebhardt 199C| , |Carlberg, Yee fc Ellingson 1997 ). As we increase the redshift 



range of the data given to the bi-weight estimator up to 25%, the velocity dispersion rises an 
average of 13% and then remains reasonably stable. We adjust the ESO velocity dispersions 
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downwards by 13%, to arrive a velocity dispersion of 708 kms -1 . We derive a mass limit for the 
ESO > 800 km s -1 sample of Mi. 5 = 4.0 x 10 14 /i _1 M . The sample is not X-ray selected so needs 
no density correction. 

An upper limit to the low redshift density for the Northern Abell sample with velocity 
dispersions has been derived previously ( White, Efstathiou &; Frenk 1993] ). The same sample, 



when compared with X-ray results confirms that the median velocity dispersion is somewhat 
overestimated ( Eke, Cole fc Frenk 1996; referred to as ECF), and argues that the velocity dispersion 



is about 650 kms -1 . Because this sample is very similar to the ESO sample, we adopt the same 
minimum velocity dispersion as we derived above, 708 kms -1 , and hence the same minimum 
mass. This sample requires no density correction. 



3.4. A High-Redshift Sample 

The EMSS sample contains a fair sample of clusters ranging from redshifts of about 0.14 to 
redshift 0.83 ( Henry et al. 1992 , Luppino Sz Gioia 1995| ). The high redshift sample ( Luppino fc 



Gioia 1995| ) is cautiously assigned the same minimum mass, M1.5 = 4.7 x 10 14 /i 1 M as we used 



for the HA subsample. Given the richness of these clusters we expect that the clusters actually 
have higher masses, consequently this is a conservative assumption. The mean density derived 
from the 4 clusters in the redshift range 0.55 < z < 0.85 for S7 = 0.2 and 1 is given in the Table 
and includes the f x correction. Reassuringly, they are similar to the densities derived elsewhere 
( Luppino &: Gioia 1995| ) using the same clusters, but a different analysis. It is known that at least 



one of these clusters, MS1054— 03, likely has a mass well over the limit to be included in the sample 
( Luppino fc Kaiser 1996j ). The minimum mass is uncertain but likely at least that of the moderate 



z CNOC sample. In particular, these are very rich clusters which the very strong richness-<7„ 



relation (Carlberg et al. 1996) indicates to be high a v clusters. We expect a v ~ 900 kms . 



4. Parameter Probabilities and Conclusions 

The x 2 probability contours of the PS predictions and the observed densities are plotted in 
Figure f|. The sample variances are calculated as the quadrature sum of the l/y/N of the sample 
and the errors in the sample correction factors. The observed densities as a function of $7 are 
estimated as a linear interpolation between the values at £7 = 0.2 and $7 = 1 measurements. The 
mass errors are known to be about 25% for a single cluster ( Carlberg et al. 1996j ), which in the 



mean for the various samples will be reduced to about 8-15%, depending on sample size. This 
relatively small error is neglected in calculating % 2 . 

The minimum \ 2 ls near $7 ~ 0.4 and as — 0.75, although Q between 0.2 and 0.6 are 
statistically acceptable within the 90% confidence interval. For our preferred $7 = 0.2 we find that 
as = 0.95 ± 0.05. A model with 17 = 1 is excluded at more than 99% confidence, although this is 
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dependent on the the estimated masses of the high redshift EMSS sample. 

We thank Mike Hudson, Jim Bartlett, Alain Blanchard and Jamila Oukbir for comments. The 
Canadian Time Assignment Committee allocated CFHT observing time. The CFHT organization 
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NSERC and NRC of Canada. 
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Table 1: Cluster Densities for 0, = 0.2 



redsnirt 


Sample 


a v (min) 
kms -1 


log M 


log 
h 


n{> M) 
" 3 Mpc 3 


0.0-0.10 


WEF 


850 


14.66 




-5.40 


0.0-0.10 


ECF 


708 


14.60 




-5.40 


0.0-0.10 


ESO 


708 


14.60 




-5.60 


0.0-0.10 


HA 


750 


14.67 




-6.12 


0.18-0.35 


CNOC 


800 


14.68 




-6.53 


0.35-0.55 


CNOC 


800 


14.80 




-7.15 


0.55-0.85 


EMSS 


800 


14.81 




-7.30 
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Fig. 1. — The n(> M) relation for the different redshift ranges as calculated from the PS formula 
for VL = 0.2 and as = 0.95. The errors are N~ 1 / 2 estimates with the selection bias errors added 
in quadrature. Optically selected clusters have open symbols, X-ray selected clusters have closed 
symbols. The CNOC and EMSS samples have densities derived in this paper. The low redshift 
samples are discussed in the text. 

Fig. 2. — The n(> M) relation from the PS formula calculated for £1 = 1 and as = 0.55. These 
model parameters do not adequately describe the data. 

Fig. 3. — The bolometric X-ray luminosity versus line-of-sight velocity dispersion for the Edge 
& Stewart and CNOC samples. The sample limits at L x > lO^ergs" 1 and a v > 800 km s^ 1 are 
shown. The ES velocity dispersions are divided by 1.13. 

Fig. 4. — A plot of x 2 for the all independent samples (solid lines) and excluding the high redshift 
EMSS sample (dotted lines). The contours are the 90% and 99% confidence levels. The results of 
the CNOC analysis, = 0.19 ± 0.06, with its la range are indicated. 
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Fig. 1. — The n(> M) relation for the different redshift ranges as calculated from the PS formula 
for 0, = 0.2 and as = 0.95. The errors are iV -1 / 2 estimates with the selection bias errors added 
in quadrature. Optically selected clusters have open symbols, X-ray selected clusters have closed 
symbols. The CNOC and EMSS samples have densities derived in this paper. The low redshift 
samples are discussed in the text. 
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Fig. 2. — The n(> M) relation from the PS formula calculated for Q = 1 and as = 0.55. These 
model parameters do not adequately describe the data. 
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Fig. 3. — The bolometric X-ray luminosity versus line-of-sight velocity dispersion for the Edge 
& Stewart and CNOC samples. The sample limits at L x > 10 45 ergs _1 and a v > 800kms _1 arc 
shown. The ES velocity dispersions are divided by 1.13. 



-13- 



C2 



00 

♦ 

o 



CD 

♦ 

O 



* 

o 



CN 

♦ 

o 



o 




0.4 



0.6 



0.8 




Fig. 4. — A plot of x 2 f° r the all independent samples (solid lines) and excluding the high redshift 
EMSS sample (dotted lines). The contours are the 90% and 99% confidence levels. The results of 
the CNOC analysis, fi = 0.19 ± 0.06, with its la range are indicated. 



